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Abstract. Preliminary results on particle spectra and fluctuations in inelastic p+p 
interactions measured by NA61/SHINE at the CERN SPS are presented and its future ion 
program is discussed. We also show results on the excitation functions of event-by-event particle 
ratio fluctuations in central Pb+Pb collisions from NA49 using a novel approach. We found 
a dependence of the results on phase space coverage which appears to explain the reported 
difference between measurements of NA49 and those of STAR in central Au+Au collisions. 
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As the fundamental theory of strong interactions, Quantum Chromodynamics (QCD), is 
asymptotically free, the created matter at high temperature and/or high density may be 
dominated by the state of quasi-free quarks and gluons referred to as the Quark-Gluon Plasma 

t-h (QGP). By colliding heavy-ions at high energies one hopes to heat and/or compress the matter 

to energy densities at which the production of the QGP begins. Indeed, lattice QCD calculations 

j — , predict this new phase at high temperatures. The motivation of investigating relativistic 

heavy ion collisions is the experimental study of hadronic matter under extreme conditions. 
A wealth of ideas in the past few decades have been proposed to explore the phase structure 
of strongly interacting matter. Among other possible probes event-by-event fluctuations of 
different observables may be sensitive to the transitions between hadronic and partonic phases. 
In particular, the location of the critical point may be signalled by a characteristic pattern in 
the energy and system size dependence of the measured fluctuation signals. 

The SPS Heavy Ion and Neutrino Experiment (NA61/SHINE) is operating since 2007 in the 

X North Area of the CERN Super Proton Synchrotron (SPS). NA61/SHINE uses the apparatus 

of the NA49 pQ experiment with numerous upgrades [2j- Its large acceptance and excellent 
particle identification capabilities [U [2] make possible a rich physics program, in particular the 
ion program which is the main subject of this contribution. The ultimate goal is to explore 
the phase diagram of strongly interacting matter using collisions of different systems from p+p, 
through p+A, to A+A collisions at projectile momenta of 13,4, 20,4, 30,4, 40,4, 80,4 and 158,4 
GeV/c. The energy scan of p+p reactions was completed in 2009-2011 and preliminary results 
will be presented. Data on Be+Be collisions at three top beam momenta and at 13,4 GeV/c were 
recorded in 2011 and 2012. In 2013 Be+Be data at the remaining energies will be collected. The 
energy scans of p+Pb, Ar+Ca and Xe+La collisions are foreseen until 2016. In addition a beam 
energy scan of Pb+Pb collisions is planned. The explored energy range probes an important 
region in the phase diagram of strongly interacting matter. Indeed, the NA49 collaboration 
reported signals for the onset of deconfinement encoded in non-monotonic behavior of excitation 
functions of several hadronic observables [3] . These observations have been recently confirmed 



by the RHIC and LHC data [I] . Moreover, the critical point at the location predicted by several 
theory groups [3] can be probed at SPS energies. 




Figure 1. (Color online) Left: Transverse- mass spectra at mid-rapidity for negatively charged 
pions fitted with an exponential function in the 0.2 < my — < 0.7 GeV/c 2 interval. Right: 
Rapidity spectra of negatively charged pions fitted with the sum of two Gaussian functions. 
Results refer to inelastic p+p interactions at beam momenta of 20-158 GeV/c. 



2. Inclusive particle spectra in p+p interactions 

Results from the energy scan of p+p interactions will serve as an important reference for 
interpretation of the p+A and A+A data at corresponding beam energies collected with the 
same apparatus. Figure [I] (left) presents the transverse mass (m-r) spectra at mid-rapidity 
for negatively charged pions in inelastic p+p collisions at 20, 31, 40, 80 and 158 GeV/c [6]. 
The results are corrected for reconstruction inefficiencies, acceptance, and feed-down using 
Monte Carlo simulations. The mj- spectra are extrapolated by fitting the distributions with 
an exponential function. Figure [l] (right) depicts the m-r integrated rapidity spectra. The latter 
allow to obtain the mean number of negatively charged pions in full phase space from which 
one can calculate the total number of pions using isospin coefficients. The NA61/SHINE pion 
multiplicities are in good agreement with the world data [6]. 



3. Event by Event fluctuations in Pb+Pb and p+p collisions. 

A motivation for our work presented in this section is to shed light on reported differences 
between the results from STAR and NA49 on fluctuations of identified hadrons [9] . The 
studied measure of dynamical particle ratio fluctuations t^^LA, £>] is defined as [lOj : 

(A(A-l)) (B(B-l)) (AB) 
Vdyn[A,B] — — 2 + — ^ 2^^, (1) 

where A and B stand for multiplicities of different particle species. 

As is evident from Eq.(l) all second moments of the multiplicity distributions of particles A 
and B have to be known in order to construct i/^yn [A,B]. The standard approach is to identify 
these particles event-by-event. However, this approach is hampered by incomplete particle 
identification, which is taken care of by either selecting suitable phase space regions or by 
applying a fitting procedure event-by-event. The latter typically introduces artificial correlations 
which are usually corrected for by the event mixing technique. 
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Figure 2. (Color Online) Energy dependence of (a) Vdyn\p+P-, 7r + +7r~], (b) Vdyn[K + +K~ ,p+p\ 
and (c) Vdyn[K + + i^~,7r + + Results from the Identity Method for central Pb+Pb data 
of NA49 are shown by blue squares. Published NA49 results are indicated by circles. Stars 
represent results of the STAR collaboration for central Au+Au collisions. 



In this context a novel approach, called Identity Method [1 lj. I12j was developed for fluctuation 
analysis. It uses a probabilistic approach which avoids the event-by-event fitting as well as has a 
rigorous mathematical derivation. Furthermore, there is no need for corrections based on event 
mixing anymore. Results using the Identity Method are presented in Fig. [2] for central Pb+Pb 
data of NA49 and compared to published NA49 and STAR results from its Beam Energy Scan 
(BES) program [9]. The energy dependence of Vd yn [p + P, ^ + + vr~] is consistent with the NA49 
publication. Moreover the increasing trend of the excitation functions of Ud yn [K + + K~ ,p + p] 
and Udyn[K + + K~,ir + + 7r~] is confirmed by this analysis. It was found in Ref. |13] that 
Vdyn exhibits an intrinsic dependence on the multiplicities of accepted particles. Published 
NA49 results are essentially performed at forward rapidities, whereas the STAR acceptance 
covers the mid-rapidity region without the low pt range. As an example we illustrate in Fig. [3] 
the phase space coverage of pions, kaons and protons at different projectile energies in the 
acceptance of NA49. In the same figure the acceptance of the STAR apparatus at corresponding 
center-of-mass energies are presented with colored lines. We, therefore, studied the acceptance 
dependence of Vdyn by performing the analysis in different acceptance bins going from forward 
rapidity to mid-rapidity. Technically different acceptance bins were selected by applying several 
upper momentum cuts for each reconstructed track and calculating the corresponding maximum 
rapidity of protons at pr=0. Thereafter we will call this quantity a proton rapidity cut. The 
dependence of obtained results on the proton rapidity cut is illustrated in Fig.|4j At 20^4 and 30 A 
GeV/c VdyniK^ + K~ ,p+p] and Udy n [K + + K~ , k + + ir~] show a strong acceptance dependence, 
and eventually hit the STAR point in a particular acceptance bin. Interestingly the acceptance 
dependence weakens above 30 A GeV/c where no difference was observed with STAR. It is also 
remarkable that Vdyn[p + p, vr + + ir~] shows little acceptance dependence. This detailed study 
of fluctuation results in different acceptance bins appears to explain the difference between the 



Figure 3. (Color Online) The phase space coverage of pions, kaons and protons in the 
acceptance of the NA49 experiment for the Pb+Pb collisions at 30^4 GeV/c (upper panel) 
and 158^4 GeV/c (lower panel) used for the event-by-event fluctuations presented in Fig. [2j The 
acceptance of the STAR apparatus at corresponding center-of-mass energies are depicted with 
colored lines. 

STAR BES and NA49 measurements. 

As mentioned in the introduction, the obtained fluctuation signals from A+A collisions 
should be compared systematically to the baseline obtained from analysis of p+p interactions at 
corresponding energies. It is however important to properly take into account trivial differences 
between A+A and p+p collisions e.g., in the size of the colliding systems. In thermodynamics the 
quantities which are proportional to the volume of the system are called extensive quantities. For 
example, the mean number of particles in a relativistic gas within the grand canonical ensemble 
is an extensive quantity. The ratio of two extensive quantities does not depend on the system 
volume and is referred to as an intensive quantity. The studied fluctuation measure Vdyn is even 
more complicated because it is inversely proportional to the system volume or to the number of 
wounded nucleons within the Wounded Nucleon Model. On the other hand the scaled variance, 
((N 2 ) — (N) 2 )/ (N), is an intensive quantity. It is therefore more intuitive to get a baseline 
from p+p interactions for the scaled variance. The scaled variance of multiplicity distributions, 
obtained via the Identity Method is presented in Fig. [5] for the p+p data of NA61/SHINE (filled 
circles) and compared to results from the EPOS model (solid lines) |14j . These are the first 
released results on fluctuations of identified hadrons in p+p interactions. Therefore, a semi- 
quantitative agreement with the predictions of the EPOS model is remarkable. More work is 
needed to interpret remaining differences between the data and the model. 

An additional complication in the experimental study of fluctuations in A+A collisions are 
unavoidable volume fluctuations from event to event. These additional sources of fluctuations 




Figure 4. (Color Online) Acceptance dependence of (a) Vd yn \p + P> n+ + (b) Vdyn[K + + 
K~,p + p] and (c) z^j/n[-^ + + A'~,7r + + ir~] in central Pb+Pb collisions of NA49 (triangles, 
squares, dots). Stars show measurements of the STAR collaboration. Results are plotted versus 
the maximum of the proton rapidity at pt=0- 




Figure 5. (Color Online) Scaled variance for (a) p + p, (b) K + + K~ and (c) ir + + ir~ 
multiplicity distributions for p+p collisions. The NA61/SHINE results are shown by filled 
circles and compared to EPOS model calculations depicted by a solid blue curve. 



may mask the fluctuations of interest. For example, the scaled variance is sensitive to the volume 
fluctuations. A set of so-called strongly intensive quantities, which within a grand canonical 
ensemble depend neither on the system volume nor on its fluctuations were proposed in Ref [JS] . 
The experimental study of these quantities is ongoing. 



4. Summary 

We presented preliminary results from NA61/SHINE on inclusive particle production in inelastic 
p+p collisions and discussed its future ion program. Event-by-event fluctuations of particle ratios 
were analyzed using a novel approach, the Identity Method. Results on excitation functions of 
Vdyn in central Pb+Pb collisions confirm an increasing trend towards lower energies, reported 
earlier by NA49. The detailed study of Vdyn reveals a strong acceptance dependence at low 
energies for the Vdyn 

\p,K] and Vd y n[^,K]- We conclude that the different energy dependence 
of Vdyn measured by NA49 and STAR ( BES program for central Au+Au collisions) is due to 
the different phase space coverage. Moreover, using the Identity Method, the scaled variance of 
the multiplicity distributions of protons, kaons and pions were reconstructed in inelastic p+p 
collisions of NA61/SHINE. These will serve as an important reference for A+A collisions. 
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